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ABSTRACT 

Aims. We aim at determining the spatial distribution of the gas and dust in star-forming regions and address their relative abundances 
in quantitative terms. We also examine the dust opacity exponent p for spatial and/or temporal variations. 

Methods. Using mapping observations of the very dense pOphAcore, we examined standard ID and non-standard 3D methods 
to analyse data of far-infrared and submillimeter (submm) continuum radiation. The resulting dust surface density distribution can 
be compared to that of the gas. The latter was derived from the analysis of accompanying molecular line emission, observed with 
Herschel from space and with APEX from the ground. As a gas tracer we used N 2 H^, which is believed to be much less sensitive 
to freeze-out than CO and its isotopologues. Radiative transfer modelling of the N 2 H^ (J=3-2) and (J=6-5) lines with their hyperfine 
structure explicitly taken into account provides solutions for the spatial distribution of the column density A(H 2 ), hence the surface 
density distribution of the gas. 

Results. The gas-to-dust mass ratio is varying across the map, with very low values in the central regions around the core SM 1. The 
global average, = 88, is not far from the canonical value of 100, however. In p Oph A, the exponent p of the power-law description for 
the dust opacity exhibits a clear dependence on time, with high values of 2 for the envelope-dominated emission in starless Class-1 
sources to low values close to 0 for the disk-dominated emission in Class III objects, p assumes intermediate values for evolutionary 
classes in between. 

Conclusions. Since p is primarily controlled by grain size, grain growth mostly occurs in circumstellar disks. The spatial segregation 
of gas and dust, seen in projection toward the core centre, probably implies that, like C^^O, also N 2 H+ is frozen onto the grains. 

Key words, interstellar medium (ISM): general - interstellar medium: individual objects: pOphA - interstellar medium: dust, 
extinction - interstellar medium: molecules - interstellar medium: abundances - Stars: formation 


1. Introduction 


It is widely accepted that the gas-to-dust mass ratio {g2d) in 
our Galaxy has a value of about 100. Of course, this is an av¬ 
erage value and it does, in particular, pertain to the diffuse in¬ 
terstellar medium (ISM). However, this value is also commonly 
adopted for dense cloud environments, such as star-forming re¬ 
gions, and the methods used to determine that parameter are de- 
scribed for instance by [Kenyon et al.| ( p~998| ), |Draine| ( |2003| ), and 
Vuong et ST] ( |2003| ). Based on observations of CO isotopologues 
and stellar extinction values in the nearby star-forming regions in 
Taurus and Ophiuchus, |Ererking et al.| ( p^82| ) established the still 


Based on observations with APEX, which is a 12 m diameter sub¬ 
millimetre telescope at 5100m altitude on Llano Chajnantor in Chile. 
The telescope is operated by Onsala Space Observatory, Max-Planck- 
Institut fur Radioastronomie (MPIfR), and European Southern Obser¬ 
vatory (ESO). 

and also based on observations with Herschel which is an ESA 
space observatory with science instruments provided by European-led 
Principal Investigator consortia and with important participation from 
NASA. 

Figures 2,11,12,16,17 and 18 are also available in electronic form 
at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or 
via http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-hA/ 


widely used molecular gas-dust relations, that is, A(^C^O) oc Ay, 
Ay < 20 mag, and ^^C^^O/H 2 ~ 8.5 x 10“^. Eor higher Ay, satu¬ 
ration makes these relations less useful. 


Young & Scoville (1991) discussed the high average g2d of 
~ 600 obtained for molecular clouds, where gas masses refer 
to HI -h H 2 . The latter was derived from CO observations and 
the dust masses were determined from IRAS observations. Eur- 
thermore, based on data for HI and HII regions, [Spitzer| ( pTTS] ) 
quoted the rather wide range of 20 to 700 for this parameter. A 
value of about 450 was found by |Liseau et al.| ( p^95] ) for L 1688 
using CS maps and 60yum IRAS data. A high g2d would also 
be consistent with the low oxygen abundance in p Oph A found 
by|Liseau & Justtanont| ( |2009] ). Eigurell of [Brinchmann et J 
( |2()13 ), which is based on a large number of galaxies, supports 
this view on larger, extragalactic scales. 

One of the reasons for these high g2d values could be that 
large portions of the dust are at temperatures consider ably lower 
than 30 K and, hence, remained undetected by IRAS ( [Young &' 
[Scoville|1991[ ). However, at temperatures below about 25 K, the 
abundance of the dominant gas tracer, CO and its isotop ologues, 
could also become highly diminished by freeze-out (e.g., [Lip|^ 
[et al.[[2013| and references therein). One might therefore select 
spectral lines of a molecule that is seen in regions where CO is 
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Fig. 1. The core p Oph A in the continuum at 350 pm, where the differ¬ 
ent evolutionary phases are identified. The position of the dominating 
radiative source S 1 is shown by the star symbol. The colour scale is in 
units of Jy beam“^ with the circles corresponding to the beam size. The 
angular scale is shown at the bottom of the figure and the orientation in 
the sky in the upper left corner. The APEX maps of this paper have their 
origin at R.A. = 16^26^27^90, Dec = -24° 23' 57'/00 (J2000.0). 


thought to be depleted. One such molecule is diazenylium, N 2 H^ 
(e.g., Caselli et al.||1999| ), which is formed from N 2 (which re¬ 
mains in the gas phase to somewhat lower temperatures than CO) 
and is removed from the g as phase mainly by CO ( [Bergin et al. 


2002 


Qberg et al. 2005). However, as pointed out by Bergin 


et al.| ( 2002| ), even N 2 H^ (and N 2 ) may become deplete d and thus 


leave very few molecular probes (perhaps only H 2 D^, [Brunken 
et al.|2014| ) to study the centres of the densest cores 


The rotational transitions (J=3-2) and (J=6-5) of N 2 H^, at 
280 GHz and 560 GHz and with upper level energies of 27 K and 
94 K, respectively, could fulfil these requirements. In particular, 
the (1-0) line at 93 GHz, having an upper level energy similar to 
that of CO, has been widely used f or dense interstellar clouds , 
including the pOphA region (e.g., |Di Francesco et al.|[2004| ). 
However, E\lk is only 4.5 K and at considerably higher temper¬ 
atures, the population of the / = 0 and J = \ levels changes only 
slowly and does as such not contribute much to the understand¬ 
ing of the overall excitation of the molecule. 


|Friesen et al.| ( |2014| ) have used the Atacama Large Mil¬ 
limeter/submillimeter Array (ALMA) to observe p Oph A in the 
N 2 H^ (4-3) lines toward a small region (< 20" x20"), but the pa¬ 
per focused on the continuum data. There are no line profiles nor 


line maps shown for these data, which, according to Friesen et al. 


( [2014| ), were either self-absorbed or resolved out by the interfer¬ 
ometer. The authors also noted that line and continuum sources 
did not coincide, neither for SM 1 nor for SM IN, a point of po¬ 
tential importance (see below. Sect. 4.4). The primary beam of 
these observations is smaller than the single-dish beams of our 
observations. 


As revealed by pure-rotational H 2 line emission, relatively 
high temperatures, 2 10^ K, are present in the interface regions 
toward the domina ting radiation source S 1 (Fig.[^. In addition, 
|Liseau et'aL \\999) found extended [O i] 63 pm, [O i] 145pm and 
[C ii] 157 pm emission, which is indicative of temperatures in ex¬ 
cess of 100 K. These phenomena probably refer to the outer lay- 
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Fig. 2. (a) Integrated intensity, f T^dv, of N 2 H+ (3-2) at 22" resolution 
with APEX and with a sampling rate of 10"/pixel is shown in the upper 
panel. The positions of the SED-class sources of Eig.[2are shown with 
their proper designations, (b) The Herschel map of the p Oph A core in 
the N 2 H+ (6-5) line that is rotated to align with the APEX data. The (6- 
5) map has a resolution of 38" and is sampled at a rate of 16"/pixel. 
Superposed, and shown as contours, is the map in the (3-2) line at de¬ 
graded resolution, (c) The orientation of the maps in the sky. 
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Fig. 3. Top: pOphA observed at 870pm with LABOCA (left) and at 350pm with SABOCA (right). Offsets are in seconds of arc and units 
are in Jy/beam. The origin, (0,0), is the same as in Fig.[^and the numbers refer to Tablej^ Middle: The SAB OCA data (left, pur ple contours), 
degraded to the 19'.'5 resolution of the LABOCA map (in colour). The colours show the IRAM map (right) of |Motte et al.| ( |1998] ), convolved to 
19'/5 resolution, and shifted to the position of the LABOCA origin. The LABOCA data are shown as purple contours. Bottom: The observed 
spectral index cr = A log /^/A log v for the wavebands 350pm-870pm (left) and 870pm-1.3 mm (right). 
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ers of the cloud. The N 2 H^ (6-5) line could be used to probe that 
temperature regime in a coherent map. 

Already early on, the actively star-forming p Oph clouds at¬ 
tracted wide attention, largely because the nearest stellar cluster 
is forming there. This is the dominant star formation mode, in 
contrast to more isolated star formation as in the Taurus-Auriga 
clouds, for example, that is contributing only little to the stellar 
initial mass function (IMF) of the Galaxy. Not surprisingly there¬ 
fore, a huge body of literature exists, and we refer to the review 
by |Wilking et al.|p008| ) and also to a few recent specific papers 
on the distribution of molecular gas and icy grains (|Shuping et al.| 
[2000t|Pontoppidan|2006t [Andre et al.|2007l [White et al.|2015| ). 

For these clouds in p Oph, the visual extinction can locally reach 
extremely high values, (max Ay >10^ mag). 

The nearby 3' (0.1 pc) p Oph A core ( [Loren et~S^[1990[ ) is 
unique in the sense that all phases of the earliest phases of stellar 
evolution are present in a very compact region (see Figs.[^an d© 
- from starless clumps. Class -1 (SM 1), via Class 0 (VLA 1623), 
Class I (GSS 30-1) and Class II (LFAM 3) to Class III (GSS 32) 


(see, e.g.. 


Comeron et al. 1993[ [Bontemps et al.[[200H [White 


[et al.|2015 ). Potentially, this fact may make it possible to follow 
the time evolution of the dust in various locations in the core, 
including grain growth and changes in composition, such as ice 
coatings onto initially “bare” grains. Such a scenario could have 
considerable effect on theories of planetary formation (e.g., the 
location of the snow line in the circumstellar disks) and help to 
relax tight timescale requirements, particularly in the outer, less 
dense regions of the disks. The status of the current understand¬ 
ing of the dust evolution in disks has recently been summarized 
by [Testi et al.[ ( |2014] ). 

In this paper, we focus on the gas and dust contents of the 
nearby star-forming p Oph A core, observed in the far-infrared 
(FIR) and sub-millimeter spectral regions (submm). For the dis¬ 
tance we adopt the accurately determined value with the Very 
Long Baseline Array (VLBA) by [Loinard et al.[ ( [2008[ viz. 
D = 120.0^42pc). The organization is as follows: in Sect. 2, our 
mapping observations of freeze-out insensitive molecules are de¬ 
scribed. T hese observations ha ve been obtained from space with 
Herschel ( Pilbratt et al. |2Q 1 0| and from the ground with APEX 
at 5100m altitude. There, maps in the submm continuum have 
also been obtained. This section also describes the reduction of 
these data, and the results are presented in Sect. 3. These are dis¬ 
cussed in Sect. 4, illustrating specifically the aspects of spectral 
line overlap and also of geometry in the radiative transfer prob¬ 
lems. We compare the results obtained with tools widely in use 
with those arrived at with our own developments. The conclu¬ 
sions follow in Sect. 5. 


2. Observations and data reduction 

2.1. SHeFl-APEX-2 and Herschel-HIFI 


2.1.1. APEX-2 

The Atacama Pathfinder Experiment (APEX) is a 12 m single¬ 
dish telescope situated at 5100m altitude in Chile. The pointing 
accuracy of the telescope has been determined at 2"(rms). We 
used the Swedish heterodyne facility instrument APEX-2 ( [Vas-[ 
silev et al.|2008[ ) for the raster-mapping observations of p Oph A 
in the J = 3 - 2 rotational transition of the N 2 H^ molecule. The 
central rest frequency is at 279515.06711 MHz and the average 
zenith optical depth was 0.04. At 279.5 GHz, the beam size is 
EWHM = 22", the main beam efficiency is T/mb = 0.73, and the 
LSB receiver temperature was 125 K. We used the East Eourier 

Article number, page 4 of[^ 


Transform Spectrometer with 73.6 kHz wide channels, resulting 
in a velocity resolution of 0.082 km s“^ per channel and a spec¬ 
tral coverage corresponding to 200kms“^. The data reduction 
involved the fitting and removal of the baselines and was made 
with the locally available software package 


2.1.2. HIFI 


The (7=6-5) line of N 2 H^ was observed at the same time as our 
Herschel H 2 O maps (2012, GT2_rliseau_l). The observational 
details regarding the partial map data of 2010 (Table[^ that were 
obtained within the KPGT_evandish_l programme were already 
described in detail by [Bjerkeli et al.[p012| ). In their Eig. 1, off¬ 
sets are relative to the (0, 0)-position at 0^2000 = 16^26”^26H, 
^2000 =-24°24'31". The more extended mapping observations 
in 2012 were made on-the-fiy, with a reference position 14' 
away. The data reduction was made within HIPE 1(Q Eor the 
final composite map, the zero-offset is at a 2 ooo = 16^26”^27?78, 
^2000 = “24° 23' 50'/3, which is slightly offset from the APEX 
origin by (-1''6, -6'/7) and which is smaller than the beam size 
at 559 GHz, that is, 38". At this frequency in HIEI-band 1, the 
main beam efficienc}0is T/mb = 0.63. 


2.2. LABOCA and SABOCA 


The continuum observations were carried out with two bolome¬ 
ter multi-pixel cameras, viz. the Large Apex Bolometer Camera 
LABOCA ( [Siringo et al.[2009[) and the Submillimetre Array, that 
is, SABOCA ( [Siringo et al.[[2010| ). The effective operating fre¬ 
quencies correspond to 870 yum and 350yum, with about 150 and 
50yum wide filters, respectivel}j^ The resolutions are 19'/5 and 
7'/5, respectively. At the two wavelengths, the maps of pOph A 
span almost 40' x40' and 8' x8', respectively (Eig.[^, and in 
both maps the centre coordinates are at R.A. = 16^26^^27?90 and 
Dec = -24° 23' 57'/0 for the epoch of J2000.0. 

The 11'4 circular field of view (EOV) of LABOCA is deter¬ 
mined by the 295 bolometers, which are arranged with 36" spac- 
ings in between them. The number of loosely spaced bolometers 
of SABOCA is 39 and its EOV is 1'5. The EOV is under-sampled 
for both cameras and is filled by repetitive continuous scanning 
of the mapped region. Eor both LABOCA and SABOCA, the 
data reductions were made with CRUSlt3 CRUSH uses the ac¬ 
tual atmospheric opacities at the time of observation. Eor the two 
data sets, average zenith optical depths were 0.42 at 870pm and 
1.07 at 350pm, with a range of 0.55 to 1.20 for the latter. 


3. Results 

3.1. Spectral line data 

The maps of main-beam corrected integrated intensity, f Lmb dv, 
of N 2 H^ (6-5) and (3-2) are shown in Eig.[^ The mapping of the 
(6-5) line with HIEI was made at an angle with respect to the 
equatorial coordinate system, so that the map had to be rotated 
to be properly aligned with the APEX data (Eig.[^). In Eigs.[^ 
and b, the shown pixel sizes correspond to the respective sam¬ 
plings, and in b, the superposed contours show the (3-2) data 
after convolution to the (6-5) resolution (38"). 

^ ftp://yggdrasil.oso.chalmers.se/pub/xs/ 

^ Herschel Interactive Processing Environment, see 
herschel.esac.esa.int/hipe/ 

3 Mueller et al. (2014): Release Note #1, HIFI-ICC-RP-2014-001, vl.l 
http: //WWW . apex-telescope. org/instruments/ 

^ http: //WWW .submm.caltech.edu/sharc/crush/download.htm 
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Table 1. Log of observations 


Instrument 

Prog/Obs ID 
ESOlHerschel 

Mode 

Frequency 

(GHz) 

Beam 

HPBW 

Observing Date 
yy-mm-dd 

fint 

(sec) 

Map Dimension 
R.A.x Dec. 

P.A. 

APEX-2^ 

090.F-9319(A) 

N 2 H+ (3 - 2) 

279.515067 

22" 

12-09-22 to 25 

65160 

4'26 X 4'32 

0° 

HIFI^ 

1342204010^ 

N 2 H+ (6 - 5) 

558.957500 

00 

10-09-01 

6552 

5'.0 X 2'0 

35° 


1342251245^ 




12-09-21 

3608 

5r8 X LO 

35° 


1342251246" 




12-09-21 

5712 

5'0 X 2'0 

125° 


1342251247^ 




12-09-21 

5712 

5'.0 X 2'0 

125° 


1342251429^ 




12-09-26 

3608 

5r8 X LO 

35° 

LABOCA^ 

090.F-9304(A) 

Continuum 

345 

19'.'5 

12-08-17 

9298 

38G X 35r5 

0° 

SABOCA^ 

090.F-9304(B) 

Continuum 

852 

7'.'5 

12-08-22 

20419 

ins X 7f85 

0° 


" Map offsets refer to 0^2000 = 16‘’26”27!*90, ^2000 = -24° 23' 57'/0. 

^ The HIFI data refer to OD475 on 10-09-01 and OD 1221 on 12-09-16 and OD 1231 on 12-09-26, with a total integration time 
of 25192 s. 

Observed centre coordinates are Q' 2 ooo= 16‘'26”26B8, ^ 2000 = -24°24'31"; 16'’26”13*79, -24°22'31"; ["] 16'’26”21?30, 

-24°26'15"; [4] 16‘’26"'3l!‘73, -24°22'47"; P] 16‘’26"'38®97, -24°26'32". 


^ 0.6 
^ 0.4 

I 0.2 
0.0 


Fig. 4. Centre of the three sub-regions of the mapped area showing the best-fit models (green) for the observed (blue) N 2 H^ (3-2) and (6-5) line 
profiles (left and right, respectively). The ra diative transfer with ALI takes the overlapping transitions of the 38 hyperfine components into account, 
with the frequencies of |Pagani et al.| ( [200^ . In the small frames, the numbers in the upper left corners are equatorial offsets in arcsec from the map 
origin (0,0). In the lower left frames, the scales are given for Tmb (K) and A\ - v - z^lsr (?^lsr=+ 3.5 km s"^), respectively (Figs. [^and|^. 
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3.2. Continuum data 


The reduced data are shown in Fig.|^ for both LAB OCA and 
SAB OCA, where intensity units are Jybeam"^ These intensity 
scales assume Gaussian beams given in Table[^ and in the figure, 
the contour levels for LABOCA are (1) -1.21, (2) 0.13, (3) 0.59, 
(4) 1.20, (5) 2.02, (6) 3.13, (7) 4.61, (8) 6.60 Jybeam-i. For 
SABOCA, the corresponding contours are (1) -1.01, (2) 0.31, (3) 
0.74, (4) 1.31, (5) 2.06, (6) 3.06, (7) 4.38, (8) 6.12, (9) 8.44, (10) 
11.50, (11) 15.55, (12) 20.92, (13) 28.03 Jybeam-i. The value 
of 1 Jy beam“^ corresponds to an intensity of about lOOMJy sr“^ 


for LABOCA and 600MJysr ^ for SABOCA, and with Icr of 
the rms noise corresponding to 5 and 30MJy sr“\ respectively. 

At the higher resolution of the SABOCA observations, indi¬ 
vidual sources can be identified, and their equatorial coordinates 
are presented in Tablej^ The rms pointing accuracy of the APEX 
telescope is about 2". The observed field at 350yum is much 
smaller than that at 870yum (Fig.|^. For quantitative applica¬ 
tions, the SABOCA data have to be degraded to the appropriate 
LABOCA values, that is, 19'.'5, and re-sampled onto a 4" grid. 
This is illustrated in the middle panels of Fig.where SABOCA 
contours are shown superposed on the LABOCA intensities. 
In addition, the IRAM 1.3 mm data of |Motte et al.| ( [T998] ) are 
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Fig. 5. Best-fit parameters from the ALI modelling of the observed N 2 H^ maps, yielding results as shown in Fig.|^ Modelled data use an intensity 
threshold of 3cr of the rms noise for both lines. Upper row: From left to right, the spatial distributions are shown for the gas column density, 
A(H 2 ), the kinetic gas temperature, Tkin, the gas volume density, w(H 2 ) and the relative molecular abundance, X(N 2 H+)=fz(N 2 H+)/w(H 2 ). Lower 
row: From left to right, the spatial distributions are shown for the line centre optical depth of the (3-2) and (6-5) lines, the N 2 H^ column density, 
A(N 2 H+), and the residuals of the radial velocity fits Av in kms"^ Spatial scales are in arcsec and the colour coding is given by the vertical bars 
next to each frame. 


Table 2. SAB OCA objects at 7'/5 resolution 


ID 

R.A. (J2000) 
hh mm ss.sss 

Dec. (J2000) 

o / // 

Compactness/Source 

Sal 

16 26 17.30 

-24 23 46.5 

diffuse, ISO-Oph 21 

Sa2 

16 26 21.53 

-24 23 04.5 

point, GSS 30-1 

Sa3 

16 26 21.75 

-24 22 51.0 

point, LFAM 1 

Sa4 

16 26 23.67 

-24 24 39.8 

point, LFAM 3 

Sa5 

16 26 26.42 

-24 24 30.4 

point, VLA 1623 

Sa6 

16 26 27.49 

-24 23 56.2 

diffuse, SM IN 

Sa7 

16 26 27.52 

-24 23 54.8 

diffuse, SM 1 

Sa8 

16 26 27.74 

-24 22 55.5 

diffuse 

Sa9 

16 26 32.68 

-24 26 09.8 

diffuse 


shown, also convolved to the 19''5 resolution of LAB OCA and 
sampled onto the 4" regular grid. For the comparison with the 
APEX results, the map of p Oph A was shifted to the same cen¬ 
tre coordinates. 


4. Discussion 

We start our discussion with determining the physical parame¬ 
ters for the gas, in particular and A^(H 2 ), using spectral line 
maps of (presumably) optically thin emission, and has of- 
ten been use d to guarantee this. However, as shown by [Liseau 


et al. 


( |2010| ) on the basis of C**0 (3-2) and (3-2) obser- 


vations of p Oph A, optical depths do exceed unity even in the 
rare isotopologues of CO. In those cases, the observed radiation 
does not carry information about the entire core but is restricted 
to layers closer to the surface, leading to erroneous estimates 
of the overall temperature and density distribution. For instance, 
temperatures based on such observations are around the (for a 



cu 



Relative radial velocity (km s'^) 


Fig. 6. Hyperfine line strengths (red vertical bars) for the / = 3 - 2 
(upper panel) and J - 6 - 5 (lower panel) transitions. The blue curves 
are theoretical profiles for an optical depth of 0.1 and a line width of 
0.1kms“^ scaled to the same maximum amplitude as the LTE line 
strengths. 


dense core) surprisingly high value of ~ 25 K. At temperatures 
below this, which are commonly encountered in dense cores, the 
CO gas will be depleted because of freeze-out of the molecules 
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onto the surfaces of the accompanying grains. In contrast, N 2 H^ 
molecules remain in the gas phase down to much lower temper¬ 
atures ( [Lippok et al.|2013| and references therein). Except for a 
few exceptions, only the ( 1 - 0 ) line and its hyperfine components 
have been observed and the analysis generally assumed ther¬ 
modynamic equilibrium (LTE). Here, we make multi-transition 
analyses and treat the radiative transfer in non-LTE. The impli¬ 
cations for our observations are described in the next section. 


4.1. Gas temperature, opacity, and mass 

Eor large numerical grids, the publicly available computer code 
radex ( [van der Tak et al.||2007| ) is very suitable. In its present 
form, the code is unable to properly handle the transfer in opti¬ 
cally thick lines with extremely closely spaced transitions, how¬ 
ever (the “line overlap” problem, Eig.|^. In addition, by con¬ 
struction, the code will not return the proper shapes of the lines. 
To surmount these difficulties, we used for the analysis an ac¬ 
celerated lambda iteration (ALI) code that takes line-overlap ex¬ 
plicitly into account and also correctly computes the line profiles 
with which the observations can be compared directly. 

The upper level energy of the N 2 H^ (3-2) transition corre¬ 
sponds to nearly 28 K, close to previous estimates of the kinetic 
gas temperature. The observed (3-2) lines exhibit the remarkable 
feature that the satellite lines reveal strengths that are almost sim¬ 
ilar to that of the main line (Eig.[^, that is, much stronger than 
expected for optically thin radiation in LTE (Eig.[^. In fact, this 
strongly suggests optically very thick emission in the overlap¬ 
ping components of the central, spectrally non-resolved, main 
line. This line anomaly is in contrast to the observed (6-5) line 
profile, which appears nearly Gaussian (however, the spectral 
resolution is lower). The J = 6 level is about 94 K above ground. 

The Onsala ALI code (P. Bergman) is based on the work by 
[Rybicki & Hummer] ( | 199 1| ) and |Rybicki & Hummer] (|1992|) and 
the bench-marking has been reported by Maercker et al.| ( |2()08] ). 
It takes line overlap explicitly into account and, in the following, 
we use this program to model the multitude of hfs lines of N 2 H^. 
The geometry is spherically symmetric, which needs to be kept 
in mind when comparing model line profiles with observed ones. 
Dust continua are included, but switched out in the simulated on- 
off observations presented below (zero baseline). 

The radiative transfer in the hfs transitions with line overlap 
was computed using the energy levels that we derived from the 
frequencies presented by [Pagani et al.| ( |2Q09| ). The collisional 
de-exc itation rate coefficients were adopted from |Lique et al.| 
( |2Q15| ), who have calculated the rates for collisions with H 2 for 
the temperature range 5 to 70 K and for the lowest hfs levels. 
However, E. Lique kindly made available to us their entire results 
for upper energy levels up to /^p = 7 prior to publication. The 
computation of the inverse rates, that is, the collisional excitation 
rate coefficients, assumed thermal balance. 

The region was devided into three parts: the north, centred on 
N 2, the middle, centred on SM 1, and the southeast, centred on 
N 6 , where the N-nomenclature refers to the source designations 
by |Di Erancesco et al. ( |2004| ). However, only pixels with > 3cr 
integrated intensities in both the N 2 H^ (3-2) and (6-5) lines were 
included. These regions, assumed to be spherical, are about 75" 
in size, hence neighbouring regions are overlapping and radii 
are 0.02 pc (6.75 x 10^^ cm). We kept the turbulent velocity con¬ 
stant throughout the core, that is, L'turb = 0.25 km s“^ After a 
number of initial trials spanning over a much larger parameter 
space, the ALI code was run on a finer grid around the most 
probable values, defined by 4.0 < log^(H 2 )< 7.0 in steps of 
0.2, -11.0 < logA(N 2 H+)< -8.0 by 0.2 and 5 < Tkin < 50K 


by increments of 1 K, generating 11776 models. Among these, 
2816 solutions were found that fit the observed intensity ra¬ 
tios within the observational errors. Using the individual line 
profiles, that is, both the shape and the absolute intensity, for 
both lines, N 2 H^ (3-2) and (6-5) from the same position (pixel), 
the best-fit ALI solution was selected for that pixel. As best- 
fit-estimator, we used sums of the residuals, where the higher 
frequency spectra were weighted more strongly to compensate 
for their overall lower intensities. In other words, we minimized 
Z3^2 (7’obs - 7’ALif + 10 X 2 : 6^5 (7’obs - ^ALif, resulting in a 
map of best-fit model profiles that, together with the observed 
ones, are displayed in Eig.|^ 

In Eig.[^ the spatial distributions of the best-fit model pa¬ 
rameters A/^(H 2 ), Tkin, ^(H 2 ) and A(N 2 H^) are shown. Values for 
the 10 " pixels span two orders of magnitude in column density, 
V(H 2 ). Eor the star pOph, [Bohlin et al^ ( |1978| ) determined the 
anomalous relation V(HI -r }l 2 )/E(B - V) = 15.4 x 10^^ atoms 
cm“^ mag“^. Since Ry = Ay/E(B - V) = 5.5, the gas-dust rela¬ 
tion for pOph reads A^(H 2 ) = 1.4 x 10^^ Aycm"^. Using this 
calibration also for the core, the ALI results over 30" scales 
would imply the surprisingly low visual extinction. Ay, through 
pOphA of about 100 mag or lower. Average column densities 
of a few times 10 ^^ cm“^ are in accord with previous determina¬ 
tions. 

Kinetic gas temperatures are about 10-13 K in the central 
regions, but 17-19 K in the north, andl6-19K in the southeast. 
Overall, these tem peratures are lower th an what has been deter¬ 
mined before (e.g., [Bergman et al.|201 1[ and references therein). 
The relative abundance, A(N 2 H^) = A^(N 2 H^)/A^(H 2 ) appears to 
peak in the central part, where logA(N 2 H^)= -9.0 + 0.2, 
whereas A(N 2 H^) is generally lower than 10 “^^ in the outer re¬ 
gions. There, derived overall gas volume densities, ^(H 2 ), are 
around 10 ^cm“^, which would be in accord with what might 
be required by some models of the dust opacity (see below. 
Sect. 4.4.2) Toward SM 1, and on much smaller scales, [Friesenj 


et al. 


(2014) estimated densities in excess of 10 ^ cm 


As expected, optical depths in the centre of the (3-2) main 
line are high, explaining the low contrast between main and 
satellite components (Eig. In contrast, the (6-5) lines are op¬ 
tically thin, which is supported by their observed Gaussian line 
shapes. The figure also shows the distribution of the radial ve¬ 
locity residuals, Av = l'o - t^LSR, where vq is the observed central 
velocity and r'LSR=+3.5 km s“^ Locally, velocity shifts are lower 
than 0.2 km s“\ which is much lower than the assumed turbulent 
velocity of the gas (0.6kms“^). On much larger scales, > 40', 
[Liseau et "SI] ( [1995 ) derived a NW-SE radial velocity gradient of 
about 1 km s ^ pc“^ . Eor the limited extent of the region studied 
here, this has only a negligible effect on the line transfer. 


4.2. Dust temperature, opacity, and mass 
4.2.1. Point-source fitting of the SED 

The fluxes of several point sources in the field were extracted 
and corrected for the underlying cloud emission. Their SEDs 
were determined assuming modified blackbody emission, where 
the adjustable parameter p is the power-law exponent of the fre¬ 
quency dependence of the flux in the Rayleigh-Jeans regime (see 
Appendix A). The fitting results are presented in Tablej^ and 
shown in Eig.|^ These suggest a temporal evolution of the pa¬ 
rameter p, in the sense that the EIR/submm SEDs become in¬ 
creasingly flatter with time. At early times. Class -1 and Class 0, 
the EIR/submm emission is entirely dominated by the envelope, 
whereas at the later stages, the envelope emission has essen- 
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Fig. 7. Left: Time evolution of the opacity exponent/5, where Ky ocv^, and the values are derived from the fits to the adjacent SEDs. Right: SEDs 
of the point sources SM 1, VLA1623, GSS 30-IRSl, LEAM1, LEAM 3 AND GSS 32, belonging to different evolutionary classes (see Eig.|^. The 
/5-values decrease from 2 for the starless dense dust core SM 1 to 0.2 for the visible TTauri star GSS 32. Colour coding of the symbols is light 
blue: ISO-CVE 5.5 pm and \6pm. Red: WISE 22pm or Spitzer 24 pm. Blue: H erschel-FACS 100//m and 160 p m White: PACS-spe ctra at 63, 73, 
145 and 1^9pm. Green: APEX 350 and ^10pm. Black refers to literature data: [Wilson et al.| ( [l999| 450//m) and jMotte et al.|(T9^ 1.3 mm). 


tially vanished and is dominated by the disk or the remnants of 
it (Class III). A mixture of these two dust-emission components 
is exhibited by the objects in between (Class I and II). The value 
of /5 is presumably dominated by the size distribution of the dust 
grains (e.g., [Miyake & Nakagawa|1993t|Kruegel & Siebenmor- 
|gen|1994] ). 

To examine the yS-dependence on the size a of (spherical) 
grains, we performed Mie-calculations ( |Min et JT][2003[ ) us¬ 
ing laboratory measurements for a few typical grai n materials 
(Fig. [8] ). This plot is qualitatively similar to Fig. 4 of [Testi et al. 


P014[ ) for young disks, but here, p has been calculated for the 
wavelength region 200/im to 1.3 mm instead, and it falls into 
the interval 0.0 < p < 2.5. The power-law exponent p has 
been varied within the range -5.0 < p < -2.0, and where 
dn{a) oc a^da. As expected, theyS-dependence on chemical com¬ 
position decreases with increasing p, that is, grain size, and these 
larger grains are found around the Class II and III objects, viz. 
LFAM 3 and GSS 32, respectively, with a 2 100yum. In contrast, 
the very small grains are found predominantly around the prob¬ 
ably youngest object, that is, the starless clump SM 1 (Class-1). 
Around the Class 0 to I sources VLA 1623, GSS 30IRS 1 and 
LFAM 1, sizes would fall roughly within 0.1 < a < lOyum. 


Therefore, the observed time dependence suggests that grain 
growth mainly occurs in the protoplanetary disks, with only lit¬ 
tle enrichment of large grains in the parental cloud. However, the 
precise dating of these events is difficult, but they probably fall 
within the range of some thousand years (starless cores. Class - 
1) to some megayears (near the ZAMS, Class III). If our obser¬ 
vational results apply to currently adopted pre-main-sequence 
timescales (e.g., [Evans et~SI][2009[ ), then pit) can be fit by a 
(broken) power law, p ~ where yi< 10^ yr) ~ 0.04 and 
yi> 10^ yr) ~ 0.5 to 0.6, respectively. If p scales inversely with 
grain size/mass, this may be consistent with stochastic processes 
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Table 3. SED-fitting result for sources at different evolutionary stages. 


Source 

7’dust [K] 

p 

Type 

SMI^ 

15 + I 

1.95+ 0.13 

Class-I 

VLAI623 

19 + 2 

1.38+ 0.14 

Class 0 

GSS30-IRSI 

65 + 4 

1.49 + 0.17 

Class I 

LFAM I 

51 + 6 

1.08+ 0.17 

Class I 

LFAM 3 

39 + 3 

0.47 + 0.21 

Class II 

GSS 32 

96+18 

>0.2 

Class III 


^ [Ward-Thompson et"ak] ( [1989| ) determined Tdust = 15 + 5 K 
andyS = 2.2 + 0.3. 


in a turbulent medium (see also [Blum][2004[ ). Apparently, very 
little happens up to some 10^ years. Thereafter, grain growth 
seems to occur at a faster rate. 

Needless to say, it would be very valuable to extend this type 
of analysis to other star-forming regions, both for better statis¬ 
tics, but above all also for the comparison among different envi¬ 
ronments. yS-evolution around young sources has recently been 
found elsewhere, but generally only for a single o bject at a time 
or, in one case, for two sources ( [Testi et al.|2Q14[ and references 
therein). However, the very compactness (0.05 pc ~ lO'^ AU) of 
an entire star-forming region in combination with the multitude 
of objects of different type as discussed here, makes the present 
study of p Oph A unique. 

As for SM 1, these results refer to source scales of some 
10^ AU. On smaller scales, that is, a fe w 10^ AU, SM 1 appear s 
to have fragmented into smaller pieces ( [Nakamura et al.|2012| ). 
This may be expected on the basis of current theoretical mod¬ 
els of star formation ( [Lomax et al. 120141 and references therein). 
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Fig. 8. Grain size exponent p, defined through dn{a) oc a^da for spheri¬ 
cal grains of size a, as a function of the grain opacity exponent p for the 
wavelength range 200 pm to 1 .3 mm, and where Ky oc yp. The curves for 
three different mat erials are shown, viz. for amorphous ic e (H 2 O, blue: 
Bertie et al.|1968| ), amorphous silicate s (MgFeSi 04 , red: [Jaeger et al.] 
1994} and aWrphous carbon (C, green: [BTanco et al.|1991 1. The values 
of p determined from observations are shown as dashed lines, referring 
to the sources of Tablej^and identified to the right. 



Fig. 9. Extinction curves used for Model I, where the mass extinction 
coefficient is given by /Cext(d) = /^abs(d) + /^scaW). Solid lines are for 
mass absorption coefficients /Cabs and dashes are for the scattering, Ksca- 
The cu rves are based on the opt i cal co nstants provided by [Jaeger et al.| 
\\994) and by jDorschner et al.] ( |1995j ) for two grain populations. The 
adopted grain size distributions have an exponent p - 3.5, and minimal 
sizes are amm = 0.03 pm for both. The largest sizes are ^max = 3 pm 
(blue) and a^ax = 300pm (red). 


However, using ALMA at higher resolution, [Friesen et al. (2014) 
found no evidence for substructures down to scales of some 
10 AU. 


4.2.2. Extended dust emission 


assume a single dust temperature nor a single power law expo¬ 
nent p of the frequency dependence of the dust opacity, nor a 
single normalization value at some fiducial frequency. Instead, 
we intend to exploit observations at, at least, three frequencies 
to determine the spatial distribution of the dust temperature, the 
dust optical depth, and the absolute grain mass opacity. However, 
from the observed spectral index a = Alog/y/Alogv for the 
wavebands 350yum-870yum and 870yum-1.3 mm, respectively, 
it is apparent that the APEX data are internally consistent, but 
appeared incompatible with the IRAM data (Fig.[^. Whereas a 
generally is > 2 for the Saboca-Laboca maps, as one would ex¬ 
pect for thermal emission in the Rayleigh-Jeans (RJ) regime, the 
index is < 2 everywhere for the Laboca-IRAM data. We have 
to conclude that the relative calibration of the LAB OCA and 
SAB OCA maps is good, but that the LAB OCA and IRAM data 
suffer from inter-calibration uncertainties, and the latter were 
therefore not used in our analysis (see Appendix A). 

With ly obtained from observation, the right-hand side of 
Eq. |A.2| contains two unknowns. For the interpretation of obser¬ 
vations of molecular clouds, it has been common to assume ad 
hoc values for Ky and In addition, one also assumes that the 
run of K with v is given by a power law, with an exponent p, that 
is, /Cy cx y For thermal emission in the RJ-regime, p would be 
related to the spectral index a, given byp = a-2, and for spher¬ 


ical Mie particles, p takes values between 0 and 2 (e.g., Emerson 


1988), but values lower than 0 and greater than 2 are possible 


when the proper physics of non-perfect crystalline material are 
included ( jMeny et al.|20Q7| ). The use of the yS-relation for mass 
esti mates requires the normalization at some fiducial frequency 
(cf. [Hildebrandj 19 8^ . In many cases, the normalization value of 
the opacity incorporates an assumed g2d of 100 ( [Beckwith et al.j 
[2000 and references therein). 


The normalization is commonly taken from published works 
on dust opacity computations. However, by using these theoret¬ 
ical results directly, one circumvents the necessity of specify¬ 
ing yS-values, which can change from one frequency band to an¬ 
other. Instead, one would wish to determine how these might 
vary spatially across a cloud, which might tell us something 
about the cloud dust properties. Potential problems include noise 
fluctuations and calibration uncertainties, and gradients in den¬ 
sity and/or temperature along the line of sight that may lead to 
non-convergence/inconsistent a:- values. Similar would apply to 
derived yS-values, when one as sumes that k{v) is given by a sin¬ 
gle power law of index p (e.g., [Juvela et al.|2Q13| ). The accuracy 
of our results below is limited to twice the statistical error of the 
intensity ly at each frequency v. However, the largest uncertainty 
comes from the absolute value of the opacity. 


As an example, we chose the work by [Qssenkopf & Hen^ 
[ning[ ( [19^ , which provides dust mass opacities for conditions 
thought to prevail in dense molecular clouds. Processes of coag¬ 
ulation are considered and the steady state is reached typically 
after some 10^ years, which is most likely younger than the ages 
of molecular clouds. The initial state consists of bare grains, dis¬ 
tributed in si ze as in the diffuse in terstellar medium (viz. “MRN- 
distribution", Mathis et al. 1977| p With time, these grains grow 


in size by acquiring ice mantles of various thickness and opaci¬ 
ties, between the wavelengths of 1 pm and 1.3 mm, are provided 
for a range in environmental density, from ^(H) = 0 to 10^ cm“^. 


The methods for estimating dust emissions from point and ex¬ 
tended sources differ only slightly (e.g., [Andre et al.|1993[ ), and 
in Appendix A, we present the analysis of our submm data as¬ 
suming conventional methods, again for the purpose of compar¬ 
ison with other works. We did not, however, as is often done, 


Draine 


^ Com paring mass a bsorption coefficients fo r MRN dust of __ 

( [2003[ ) wit h those of [Qssenkopf & Henning] \\994) indicates some 
differences: [Draine[ ([2003b lists values of = 1.92cm^g~^ and 

rh( 




nmg 


The corresponding data by Qssenkopf & Hen- 
19941) are > 3.64 cm^ g"^ ^ n 70 rr-i 


0.383 cm^ g 


’ g ^ respectively. 
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Table 4. Dust masses 


for opacites from Ossenkopf & Henning ( |1994[ ). 


Table 5. Stellar parameters used in model I 


K 

curve 

ice 

coating 

«(H) 

(cm 

^theo 

350-870pm 

Average 

rd(K) 

lOOMd 

(Me) 

1 

no ice 

0 

1.84 

21.06 

15.92 

2 


10 ^ 

1.50 

45.88 

3.28 

3 


10 ® 

1.28 

81.44 

0.81 

4 


10 ’ 

1.12 

107.06 

0.23 

5 


10 ^ 

1.05 

112.52 

0.14 

6 

thin 

0 

1.89 

19.38 

10.54 

7 


10 ® 

1.86 

20.16 

7.42 

8 


10 ® 

1.86 

20.10 

5.79 

9 


10 ’ 

1.86 

20.00 

4.97 

19 


10 ^ 

1.86 

20.24 

4.66 

20 

thick 

0 

1.89 

20.36 

7.69 

21 


10 ® 

1.89 

19.36 

6.17 

22 


10 ® 

1.89 

19.43 

5.52 

23 


10 ’ 

1.89 

19.44 

5.41 

24 


10 * 

1.92 

18.91 

5.61 


Shown below, in Sect. 4.2.4 (Fig.p^, are the results for fix¬ 
ing the opacity curve, that is. No. 22 in Table|^ The distribution 
of the dust temperature exhibits an average value of 19 K. Al¬ 
though temperatures are around 20 K in most parts of the map, 
regions with Td > 25 K exist along a north-south ridge, in ad¬ 
dition to a couple of localized shocked regions, near VLA 1623 
and GSS 30 (cf. Fig.[^. The dust optical depth at 350 pm, 
is much lower than unity, justifying a posteriori the assumption 
of low optical depth. This is then of course also simultaneously 
justified for the LABOCA data (maxrgvoyum = 0.005). 

A major drawback of the presented analysis is that it ig¬ 
nores the details of the processes that are heating the observed 
dust. In other words, the energy balance equation is not solved, 
which neglects the existence of gradients in density and radiation 
field strength, and hence temperature, in the source. These short¬ 
comings are considered in the physical models described below, 
which are two 3D radiative transfer models that follow two dif¬ 
ferent approaches. Regarding spatial coverage, these models are 
complementary, meaning that model I provides a larger, global 
view and model II focuses on the central regions of p Oph A. 


4.2.3. 3D dust modelling: Model I 

Model I uses the code RAD MC-3Cp] deve l oped b y C.R Dulle- 
mond and benchmarked by |Pinte et ah] ( |2009| ). The model 
consists of a low-density widespread cloud medium, the high- 
density core p Oph A with the dust clump SM 1, the two point 
sources VLA 1623 and GSS 30, and the two early-type stars S 1 
(B4) and HD 147889 (B2.5), each with its own PDR bubble. 
These are the main energy sources in the model, providing the 
heating of the dust. The stellar SEDs are based on Kurucz AT¬ 
LAS 9 model atmospheres with the parameters of Tablej^ 

In detail, the model assumes a high-density core of radius 
0.045 pc that is embedded in a low-density cloud. The early-type 
stars S 1 and HD 147889 generate two PDRs, viz. west and east 
of the dense core, with respective radii of 0.06 pc and 0.5 pc. 
The bright B2 star HD 147889 is situated more th an half a par¬ 
sec aw ay to the southwest and behind the cloud ( [Liseau et al.| 
1999| ). Together, these stellar radiation sources are heating the 


^ http: //WWW .ita.uni-heidelberg.de/dullemond/software/ 
radinc-3d/ 
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Parameter 

SI 

HD 147889 

M(Mo) 

6 

12 

R(Re) 

3 

7 

TeS (K) 

17 000 

22000 

log g (in dyn cm ’) 

4.0 

4.0 

I'turb (kms"') 

2.0 

2.0 


dust of p Oph A from the outside. By comparison, internal heat¬ 
ing sources are essentially negligible. 

The exploited grain opacities are shown in Fig.[^ and in 
pOphA, individual sources contribute to the total modelled 
dust mass of 10.4 x 10“^ M© as follows: the core with 7, 
SMI with 1.7, VLA 1623 with 0.9, and GSS30 with 0.8 in 
units of 10“^ Mq. a match to the PDR around HD 147889 as¬ 
sumes 0.04 Mq of dust (see Fig. 10). Based on and other 
molecules (Paper II, in preparation), the density profile is given 


by p cx for pOph A globally and by p cx r~^ for SM 1 
locally. 

In the model, the stellar object S 1 ( Chini|l9^ B5|^is sit¬ 
uated off-centre in a low-density bubble of radius 0.06 pc that 
is filled with 7 x 10“"^ M© of 0.1pm sized carbon particle^ 
( [Dorschner et al.||1995| ). These small particles scatter the short- 
wavelength radiation extremely well, and this is needed to ac¬ 
count for the observed intensity distribution in the near- and mid- 
IR. However, some finite, non-zero filling of the bubble is also 
required for numerical reasons. For an assumed g2d = 100, the 
total mass from model I of the entire pOph A core is 23.7 Mq 
(Fig.[g . 

As already mentioned, we used RADMC-3D to solve for 
the temperature distribution, which is shown along two differ¬ 
ent cuts on the right-hand side of Fig.[^ Obviously, the grain 
temperature also depends on its size, and the distributions for 
the three groups of grains are shown individually. For pm- and 
mm-size grains, temperatures of the dust amount to 10-30 K near 
the dense core SM 1. These model results can be used to gen¬ 


erate synthetic images of pOphA. These are shown in Fig. 12 


together with the observations at six wavelengths, from the mid- 
infrared to the submm. Evidently, model I reproduces the ob¬ 
served intensity distributions very well, in particular the large- 
scale features. 


4.2.4. 3D dust modelling: Model II. 

Thermal equilibrium models 


The modelling follows the procedure outlined by 
et al.| ( |20lT] ), which starts by deducing the 350pm optica 
distribution for a fixed = 2 from the observed APEX (350pm 
and 870pm) maps. The optical depth map is used to set the 
column density distribution of the model. However, the line 
of sight (LOS) density distribution is not directly constrained 
by the observations. We assumed a Gaussian distribution with 


Juvela 

depth 


This is a bina ry (B4+K, Gagne et al. 2004| . The magneti c B star ([An¬ 
dre et al.|1991|) is also a brigh t source of X-ray emission ( IHamaguchi 


et al.|2003 Gagne et al.[ 

^ For MIL - 4.6 X 10“"^ of an equivalent single star, the blow-out 
grain size is larger than 50 pm, and radiation pressure would empty the 
12 000 AU region of carbonaceous particles of size 0.1 pm within 3 500 
years. The likely existence of an intense stellar wind, exerting additional 
outward momentum, would shorten this timescale even further. 
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Fig. 10. Left panels: Density distribution of model!. With the logarithmic colour bars on top, a 2D distribution is shown to the left that lies in 
a plane containing the star S 1 and the centre of the main core SM 1, and to the right the entire model, also containing HD 1478 89. The panels 
below show the cuts along the straight white lines. White contours refer to the (3-2) measurements by Liseau et al. |( |2010| ). Right panel: 
Resulting temperature distribution of model I, with similar format as for the density distribution. The different colours in the frames below refer to 
the temperatures of the various dust constituents - Red: Very small grains (VSGs) and PAHs. Green: grains with sizes between 0.03 - 3 pm. Blue: 
grains whose sizes extend up to 300//m. 


Table 6. Assumed submm/mm dust opacities, normalized to g2d - 100. 


Opacity reference 

^350yum 

(cm^g“5 

^870yum 

(crn^g-') 

^1.3 mm 

(cm^g-i) 

(Mo) 

(kM)I(koMo) 

350, 870, 1300jum 

Note 

ID-models 

This work (Mo, 

'fo) 

11.78 

2.148 

0.9621 

5.5 

1 .0, 1.0, 1.0 

r=10-50K, OHNo.22,yS= 1.89 

Johnstone et al. 

20001 


lQa,b 


8.1 (5.2) 

0.7 (0.4) 

r=const=20K (r=13 - 42 K) 

Andre et al. (|20( 

2 ) 



0.5" 

7.5 

0.5 

T =const=12 or =20 K 

3D-models 

Model 1^ 

2.87 

0.406 

0.1846 

10.4 

0.5, 0.4, 0.4 

MgFeSi grains for = 3 pm 

Model 11.1 ^ 


11.17 

1.825 

0.8206 

7.3 

1.3, 1.1, 1.1 

Modified OH No. 22, i.e.,yS = 1.99 

Model 11.2" 


2.05 

0.376 

0.1886 

36.7 

1.2, 1.2, 1.3 

Mass for g2d = 100 

Model 11.3 ^ 


2.00 

0.393 

0.2073 

36.2 

1.1, 1.2, 1.4 

Ry = 5.5 and PAHs 


^ For the p Oph A core with assumed gas-to-dust mass ratio of g2d = 100. 
^ Refers to 850yum. 

^ Opacities from [Jaeger et al.| ( [l994| ); [Dorschner et al.| ( |1995 ). 


^ Modified Curve No. 22 of Ossenkopf & Henning (19941 using submm-yS = 1.99, but normalized at 100 pm, i.eySsso-ioOyum = 2.18. 
^ M = 59.1 Mq\ unmodified BARE-GR-S assumes g2d = 161 ( [Zubko et al.|2004| ). 

f M = 38.0 Mq\ Ry = 5.5 curve of Weingartner & Draine ( 200l| that includes PAHs assumes g2d = 105.1. 


a FWHM of 0.05 pc towards the densest part of the core, 
increasing up to 0.07 pc towards the low-density regions. 
With these values the LOS extent of the core is similar to its 
plane-of-sky dimensions. The resulting 3D density distribution 
was further modified to include a cavity around the S 1 star. 
The cavity was modelled as a paraboloid with the star S 1 at the 
focus. Inside the cavity the density was set to 10“^ times the 
highest density of the model core. The models represent a box 
with a size of 0.225 pc at a distance of 120 pc and use 256^ 


cells, yielding a resolution of approximately 180 AU, or 1'.'5. 
The three-dimensional density distribution and the plane-of-sky 
surface density distribution are illustrated in the bottom left and 
right panels of Fig.[^ respectively. 

The model was embedded in an interstellar radiation field 
with intensity three times that of the solar neighbourhood. How¬ 
ever, irradiation of the core was dominated by the individ¬ 
ual stars, mainly S 1 and HD 147889. The stellar object S 1 
was modelled as a 16 000 K, L = 1 lOOL© blackbody, while 
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Fig. 11. Images of the full 3D dust model I at 22 pm and 160 pm for the entire p Oph A core. Upper panel: Synthetic images at 22 pm and 160 
pm. Lower panel: The observations with WISE (22 pm, left) and PACS (160 pm, right). 
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HD 147889 was given a temperature of 20 000 K and luminos¬ 
ity L = 4 500 L©. Because HD 147889 is more than 0.5 pc away 
from the dense parts of the core compared to less than 0.1 pc 
for SI, the latter is the main heating source in most parts of 
the model core. In addition, the model included two deeply em¬ 
bedded stellar objects with luminosities of approximately 2Lq, 
but these do not have any significant effect on the radiation field 
except in their immediate surroundings. The radiative transfer 
calculations were made using the CRT program ( [Lunttila & Ju-| 
|vela|2012] ). The program calculates the dust temperature taking 
into account the dust self-absorption and heating, and integrates 
along the line of sight to produce maps of dust FIR emission. 
For comparison with observations, the synthetic maps were con¬ 
volved with Gaussian kernels corresponding to the APEX beam 
sizes. 


We have calculated grids of models with varying dust prop¬ 
erties and density scaling, but keeping the same radiation sources 
and the three-dimensional density distribution. In model II. 1 
we adopted the mass extinction from curve No. 22 in Table|^ 
taken from [Qssenkopf & Henning ( [1994 ); however, the opac¬ 
ities are given only for wavelengths X > 1.0 pm and no in¬ 
formation is given about scattering properties. Therefore, for 
the shorter wavelength extinction, the model Milky Way A for 
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Ry = 5.5 of Weingartner & Draine ( 200I| ) was used, matched 
to the former at the wavelength of 1.0 pm. The scattering prop¬ 
erties were obtained by assuming the same albedo and scatter¬ 
ing asymmetry parameter as in the [Weingartner & Draine| ( [2Q0r 


Ry = 5.5 model. To study the effect of jS, in addition to the un¬ 
modified p ^ 1.9 dust, we used models where the extinction 
cur ve at wavelengths A > 100//m was modified by multiplying 
the [Ossenkopf & Henning] \199A\ absorption cross-sections by 
{AllOOpmY with a = -0.1, 0.1, 0.2, and 0.3, corresponding ap¬ 
proximately lo p = 1.8, 2.0, 2.1, and 2.2. The extinction curves 
are shown in Fig.p^ and b. The peak dust surface density was 
varied in eight steps between 0.00501 g cm“^ and 0.0130 g cm“^. 


In model II. 1, the best fit was obtained with p = 2.0 and 
a peak dust surface density Edust = 0.010 gcm“^. Dust models 
with lower p predict a too high 870 pm to 350 pm brightness ra¬ 
tio. Assuming a gas-to-dust mass ratio of 100, the mass of the 
main part of the core (defined as having at least 5% of the peak 
surface density and corresponding roughly to the extent that was 
used for the ID analysis in Sect. 4.2.2) is 7.3 M©. The mass in¬ 
cluded in the whole simulated region is approximately 20 per¬ 
cent higher. Although the best models fit the observed 350yum 
and 870yum intensity within approximately 20 percent, none of 
the models is successful in reproducing the lOOyum and 160yum 
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Fig. 12. Top panels: Synthetic images of pOph A that were generated from the results of model!. Continuum wavelengths are 22, 63, 100, 160, 
350, and 870pm. The logarithmic scales are in units of Jy sr"^ and given along the colour bars atop each frame. Bottom panels: The displayed 
observations are taken from WISE (22 pm), Herschel-?AC^ (63, 100 and 160 pm), and APEX (350 and 870 pm). 
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Fig. 13. a. Modification of the extinction curves of |Ossenkopf & Henning] |I9^ (OH94) for various gr ain sizes and yS-values (see the text). The 
extinction coefficients Afabs and /Csca of model Milky Way A for Ry - 5.5 of |Weingartner & Draine| ( |2001| ) match at 1 fim. b. The extinetion curves 
of |Zubko et al.| ( [2004| ZDA04 in red) and |Weingartner & Draine| ( |2001| WDOl in black). The modified curve of |Ossenkopf & Henning | ( [T994| ) 
is shown in blue (OH94 mod). Solid lines depict mass absorption, dashed lines scattering, c. Structure of model II, where the isodensity surface 
corresponds to 20 percent of the highest density. The surface encloses approximately 30 percent of the total mass of the core. The line of sight 
(LOS) is along the negative z-axis, i.e. the axis pointing towards the observer. In the sky, right ascension is along the negative v-axis and declination 
along the positive y-axis. The dominating heating source S 1 is situated near the fo cus of the paraboloid low-den sity region, d. The surface density 


distribution of the dust, E^ust in g cm for model II.3, i.e. for the extinction curve (Weingartner & Draine 


2001 


Ry = 5.5 with PAHs). 


Herschel-FACS observations. The calculated surface brightness 
is at least a factor of two higher than observed. In particular, it is 
difficult to simultaneously produce the relatively high brightness 
at 350yum and the low observed intensity in the PACS bands. 
Figure[^ shows the comparison between the simulated maps 
from the best-fit model II. 1 and the Herschel-FACS lOOyum and 
160yum, and the SABOCA 350yum, and LABOCA 870yum maps. 
To study the effects of numerical resolution on the results, we in¬ 
creased the resolution by a factor of four, that is, to a cell size of 
0.0002 pc = 45 AU. The differences between lower and higher 
resolution were lower than 10% everywhere and lower than 2% 
where the processes were active. In the following, we therefore 
ran the models with the lower resolution to keep the computa¬ 
tions tractable. 
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In model II.2 we have chosen BARE-GR-S of IZubko et al.l 
(2004 Table 7) as our basic dust opacitie^ As with model II. 1, 
we ran simulations with modified dust models with different val¬ 
ues of p. The modifications of the extinction law were similar 
to model II. 1, except that the changes were applied at wave¬ 
lengths T > 20yum. In addition to the standard p ^ 2.0 model 
we used p = 1.5, 1.8, 2.1, and 2.2. The density scaling was var¬ 
ied with nine steps corresponding to peak dust surface densities 
(at E/5 scale) betw een 0.0197 gcm“^and 0.0985 gcm“^. As the 
[Zubko et ar] ( |2004| ) dust models have much lower sub-mm opac¬ 
ities than the dust used in model II. 1, in both cases the models 
span a similar range of 350yum optical depths. The best-fit II.2 


These opacities do not require any ad h oc adjustment s of the silicon 
or other abundances (see the discussion by Draine|2009] ). 
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Fig. 14. Results from model 11.1, comparing modeling results (left) with the observations (right). From top to bottom PACS 100 pm, PACS 160 pm, 
SABOCA350yum, and LABOCA 870//m, where intensity units are Jy sr"^ Including the PAH opacities did not improve the models at the FIR 
wavelengths, where the models over-predict the intensities observed by Herschel. In contrast, the agreement between models and APEX observa¬ 
tions in the submm is very satisfactory, within better than 20%. 
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model has a peak dust surface density Sdust = 0.0498 g cm“^ and 
P = 2.1. Modified dust with a steeper sub-mm extinction law 
matches the 870/im to 350yum brightness ratio better than the 
original dust model in this case as well. Using t he same Ay- 
to-gas mass ratio as in the BARE-GR-S model of [Zubko et al. 
P004| ) for our modified p = 2.1 dust, the total mass of the best- 
fit model is 59.1 M©, where a g2d = 161 has been assumed. 


In all models, the amount of dust in the cavity is very low, 
resulting in a visual extinction between the exciting source 
S 1 and the inner edge of the cavity of Ay = 3 x 10“^ mag. 
The extinction toward the observer is less well constrained, 
but model II.2 gives Ay 5 15 mag, with similar values for the 
others. This value compares very well with the Ay = 12.5 mag 
determined by jWard-Thompson et al.|(r989| ) from the H 2 0-ice 
absorption feature at 3.1 yum, and the range Ay = 11 to 13 mag , 
for Ry = 5.5, based on the X-ray analysis by Vuong et al. ( 2003| ). 


Non-equilibrium models 


Models II. 1 and II.2 assume that the dust grains are in 
thermal equilibrium with the local radiation field. However, 
absorption of a single high-energy photon can heat a very small 
dust grain to a high temperature. Absorbed energy is therefore 
reradiated more in the infrared than in the sub-mm (see, e.g., 
|Draine|2()03] ). To examine whether an accurate calculation of the 
emission from very small grains changes the results, model II.3 
uses the dust model Milky Way A for Ry = 5.5 of jWeingartnerj 
& Draine] ( [2001 j ) but does not assume thermal equilibrium. Be¬ 


cause these calculations are much more time-consuming, only 
three different density scalings were used with this dust model. 
At wavelengths A > 200yum the calculated emission from the 
best-fit model is almost identical to the best-fit model II.2. At 
lOOyum and 160yum, model II.3 produces slightly lower fluxes 
as shown in Fig.[T^ but the discrepancy between the model and 
Herschel-VACS observations remains large. The best-fit model 
has a peak dust surface density of Edust = 0.0491 gcm“^ and, 
using a gas-to-dust mass ratio of 105.1, a mass of 38.0 M©. 

The line-of-sight mass-averaged dust temperature, Tdust, and 
the 350yum optical depth, from the best-fit ID model (see 

Sect. 4.2.2) and the three 3D models II. 1, II.2, and II.3 are shown 
in Fig.[^ The results from all three 3D models are similar. Their 
LOS optical depths agree within 20 percent and the calculated 
dust temperatures only differ by approximately 1 to 2 K. The 
large differences in calculated masses are caused mostly by the 
widely varying sub-mm opacities and, to a smaller degree, by 
the gas-to-dust mass ratios of the dust models. The LOS aver¬ 
age temperature through the dense core is between 18 and 20 
K, but along the borders of the low-density cavity around the 
S 1 star the temperature rises to 25-30 K. Results from the ID 
model show much larger variations in temperature, although the 
average value is similar. The 350yum optical depths from the ID 
analysis are approximately 30 % lower. 


4.3. The g2d in p Oph A 

Tablej^ compares the results of different approaches to the dust 
modelling. Derived absolute masses can vary quite substan¬ 
tially, from 5 to 40 M©. However, opacity-weighted, normalized 
masses agree to within 10% of each other, with an average value 
of 1.1 over an interval of 0.4 to 1.4. 

To determine the absolute gas-to-dust mass ratio, g2d, is 
therefore not a trivial task. For example, if we choose the re¬ 
sult of the non-equilibrium model 2.III, the map average is only 
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Fig. 15. SEDs of yoOphA. For models II. 1 (black) and II.3 (blue), the 
emerging specific intensity has been integrated over the core p Oph A 
that is displayed in Fig.p^d. The t hermal equilibrium model II. 1 uses 
the extinction curve as modified by jOssenkopf & Henning |(|1 994 j). The 
non-equilibrium model II.3 uses that of jWeingartner & Draine| ( |2001j ) 
for Ry - 5.5, which includes the PAHs, and where by comparison, the 
100 and 160yum fluxes are reduced. 


slightly sub-canonical, g2d = 88 (Fig.p^. This value would di- 
rectly agree with that derived from X-ray work by j Vuong et ^ 
( 2003 ), who found g2d = 95-80 within the range Ry = 4 - 6.[^ 

However, on scales of the individual sub-regions, this num¬ 
ber would not be truly representative: Over tens of seconds of 
arc, g2d takes median values of 145 ^|q 5 around N 2 in the north, 
of 14'^Iq in the central regions of SMI, and of lOd^^^ in the 
southeast, near N 6. Taken literally, local g2d values in p Oph A 
imply fractions of the dust mass from lower than 0.5% to more 
than 10%. In particular, the deviating, low values around the core 
SM 1 were already remarked upon by [Bergman et al.| ( [2011| ). 

To bring the central g2d into line with the widely accepted 
value of 100 would require the downward adjustment of the de¬ 
rived gas mass by a factor of ten. This is not a realistic option 
on the basis of the ALI modelling, however, and must be dis¬ 
carded. Similarly, reducing the dust mass by an order of mag¬ 
nitude would result in g2d of several thousand in the north and 
southeast, an option that is not convincing either. On the other 
hand, an uneven projected distribution of gas and dust may be 
a solution. In Fig.[^ maps of gas and dust at 7" resolution are 
displayed, from which it becomes apparent that the continuum 
radiation peaks at positions offset from those of the gas. A low 
g2d would naturally result. 


4.4. Distribution of gas and dust in p Oph A 

The C^^O (3-2) line and the LAB OCA continuum band are close 
in frequency, at about 850 pm wavelength, and the maps of both 
have been observed with the same telescope (APEX), meaning 
that these data share the same angular resolution of ~ 20". Our 
SAB OCA continuum map at 350 pm has also been obtained with 
APEX, albeit at the higher resolution of 7". Both APEX contin¬ 
uum maps share the same offset centre and internally agree ex- 
celle ntly well (Sect. 4.2 ). The nominal offset of the C^^O (3-2) 
map ( [Liseau et al.|2010[ ) is at (-45'/1, -r3'/0), which was applied 


Other parameters include the metallicity Z/Z© = 0.9 - 1.0 and the 
minimum grain size amm = 0.025 - 0.07 pm. 
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Fig. 16. Distrib utions of dust temperature (left) and 350yL/m optical depth (right). Top: ID model using opacity curve No. 22 of |Ossenkopf &| 
|Henning| ( |1994| TableB.l), where p - 1.89. Second from above: 3D model (cf. Fig.[^ accounting for gradients in the core and using a slightly 
modified curve No. 22 with p - 1.99 (model II.l). The contours corres pond to 16, 17, 18, 19, 20, 25, and 30K. Second from below: The same 
3D model geometry, but with the opacity curve of [Zubko et al.| ( [2004]| that includes PAHs and f or p - 2.1 (model II.2, see the text). Bottom: 
Model II.3 has the same geometry as before, but uses opacities from jWeingartner & Draine| ( |2001| ) for Ry - 5.5 and PAHs (see Fig.p^. 
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to the line data to align with the continuum maps. As already 
mentioned, the pointing accuracy of the APEX telescope is 2" 
rms. 


We applied the program IMAGINE, which is based on an 
in-house developed algorithm ( |Rydbeck||2008| ), to the spatially 
oversampled C^^O(3-2) data. The number of noise channels in 
the line spectra, in each positional pixel, is large, which allows 
an accurate description of the statistics of the noise distribution. 
In addition, the spectra have a high S/N. These properties of the 
data are used by IMAGINE to deconvolve high-quality data, in¬ 
creasing the resolution by up to a factor of three. The decon¬ 


volved (3-2) map of the integrated line intensity f Tp^dv for 

the gas has then effectively the same resolution as the 350//m 
continuum data for the dust (Eig.[^. It is gratifying to notice 
that the point source VLA 1623, which is detected in both data 
sets, is at the same position in both maps. Eurthermore, as a fi¬ 
nal sanity check, we degraded the high-resolution result with the 
APEX beam at 850 yum, and this procedure did indeed recover 
the observed (input) map. We are therefore confident that any 
differences in the intensity distributions and on spatial scales ex¬ 
ceeding 3(T (= 6") of the pointing accuracy are real. 


This gas-dust dichotomy is not immediately apparent in 
the presented maps involving N 2 H^, since these have a much 
lower resolution, governed by the (J=6-5) observations with the 
roughly 40"-beam of Herschel. However, the (3-2) data have 
been obtained with APEX and for the same coordinates as the 
LA-/SABOCA maps. Using IMAGINE for the N 2 H+ (3-2) map 
then yields a similar resolution as the SAB OCA image and the 
sharpened (3-2) map. The result is shown in the right panel 
of Eig.[^ The clear offset, by 12", between the continuum 
source VLA 1623 and an intensity peak in N 2 H^is at first as¬ 
tonishing . Similar mismatches are also found elsewhere in these 
images, and they have also been noticed by others ( |Di Erancesco| 
|et aL|2004||Eriesen et al.|2014| ). We therefore conclude that the 
spatial non-coincidence between moleculer gas and dust is real. 
However, since the N 2 H^ emission ridge is much narrower than 
that in C^^O, a depression in the distribution is not evidenced by 
the present image. This may require data at higher resolution. 

Erom Eig. it is immediately evident that the dust con¬ 
tinuum and the molecular line emission have different spatial 
distributions in the sky. Maximum line emission in C^^O(3-2) 
is observed in a ridge closest to the stellar heating source S 1, 
whereas the continuum is mostly pronounced offset to the west, 
where all dust models, both ID and 3D, indicate dust tempera¬ 
tures Tdust ^ 20 K. In contrast, the spectral line analysis points 
toward gas temperatures that are clearly much lower than twenty 
Kelvin (see also, e.g., |Stamatellos et al.|2007| ). This might effec¬ 
tively have created the “hole” of the emission, which might 
be due to depletion of the CO gas onto the dust. At r^in < 25 X, 
carbon monoxide is mos tly frozen out ( [Sandford & Allamandola| 


1993 


Qberg et al. 2005), whereas diazenylium is possibly much 


less affected because it mainly remains in the gas phase. It is 
not clear, ho wever, how N 2 H^ w ould react at temperatures even 
below 15 K ( Qberg et al.|2005t . A selective freeze-out scenario 
would not explain the low g2d that is based on N 2 H^ observa¬ 
tions. On the other hand, if all molecules are frozen out, there 
would be no contradiction (3 


Sandford et al. { |l993t found evidence for frozen H 2 in the 
p Oph cloud. However, this result was retracted in 2000, Sci., 287, 976. 
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5. Conclusions 

In the EIR/submm, observed point sources exhibit a dependence 
of the dust parameter p on evolutionary status (SED class), 
p = p{t). Values decrease monotonically with age. Eor the pre- 
stellar core SM 1, the continuum emission is entirely dominated 
by the thick envelope and P = 2, whereas in pre-main-sequence 
stellar objects of Class III, the envelope is lost and the emission 
originates only from the circumstellar disk, yS ~ 0. Objects of 
intermediate type display intermediate yS-values. If not merely a 
matter of contrast, this may imply that grain growth primarily 
occurs in the protoplanetary disks, not in the natal cloud cores, 
since low values of p are ascribed to large grains. 

pOphA is exposed to uneven irradiation by stellar sources, 
heating the core primarily from the outside. We examined differ¬ 
ent dust emission models, including both ID and 3D geometries, 
without and with thermal balance and with non-equilibrium con¬ 
ditions. Mass estimates span a factor of seven, but do agree to 
within a factor of two when proper weighting by the absolute 
mass opacity Ky is applied. 

To model the observed gas emission in N 2 H^ lines, best-fit 
solutions were obtained by combing a multi-dimensional param¬ 
eter space. These provided excellent fits to the spatial distribu¬ 
tion of the line intensity and the line shape. Combining the re¬ 
sults for the gas with those for the dust yielded a global gas-to- 
dust mass ratio that is only slightly sub-canonical, g2d = 88. 
However, on local scales, this hardly applies anywhere, with ra¬ 
tios much higher than 10^ in the north and south and signifi¬ 
cantly lower than that in the central regions. There, gas and dust 
appear to be spatially segregated. At high angular resolution, a 
hole of C^^O (3-2) emission is accompanied by maximum emis¬ 
sion in the dust continuum. At lower resolution, this may result 
in a misleading g2d, and perhaps indicate that N 2 H^ is frozen 
out as well. 
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next to each image. 
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Appendix A: Estimation of dust masses in 
homogenous and isothermai media 

The modification we made to the commonly adopted procedure 
is based on the fact that we do not a priori assume only a single 
dust temperature and grain opacity. Rather, these quantities are 
parameters to be fit iteratively by comparison with submillimeter 
observations. The value of the optical depth at any frequency 
is obtained and tested a posteriori. The specific intensity /y, in 
erg cm“^ s“^ Hz“^ sr“\ along the line of sight of optically thin 
thermal emission from dust particles can be written as 


where Ky is the dust mass opacity, in cm^ g“^ at the frequency v, 
in Hz, is the volume density of dust particles, in cm“^, is 
the mass of a dust grain, in g, and By{T^) is the Planck function, 
in erg cm“^ s“^ Hz“^ sr“^ at v and for the dust temperature Td, 
in Kelvin. 

For constant Ky, and along a ray, the integral is only 
over the volume density meaning that it is the column density 
of dust grains, ^n^ds = N^. The surface density of dust particles 
is Ed = md A^d- Identification in Eq.jA.l [yields 




ly 

Ky ByiT^) ■ 


(A.2) 


From Eq. |A.2| it is apparent that the optical depth through the 
dust is given by 


Ty — Ky Ed- 


(A.3) 


The assumption that tv 1 will have to be justified a 
posteriori. Applying Eq. |A.2| to observations lyi at a given fre¬ 
quency vi, together with Kyi from a particular opacity curve 
(7 = 1, 2, 3, ....), yields a multitude (/ = 1, 2, 3, ...) of solutions 
for both Ed and that is. 


Sd(L7) 


lyl 

Ky\{j) Byi\T^{Uj)\' 


(A.4) 


This multitude is limited by the range of examined temper¬ 
atures, viz. Td = 4 to 2000 K, where the first value refers to 
the average equilibrium temperature of grains irradiated by the 
interstellar radiation field in the solar neighbourhooc|^ and the 
second to approximate dust evaporation temperatures. Similarly, 
the number of discretized temperature values, rd(/, 7 ), that need 


ly = ^ Ky Tid md By(T^) ds. 
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For p Oph A, the ( anisotropic) radiati on field is more intense by two 
(A. i) orders of magnitude ( [Liseau et al.|1999| ). 
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to be computed is ultimately limited by the measuring accuracy 
oily. 

This multitude can be further reduced, in fact to unity, since, 
for any given opacity cur ve, on ly one value of Td (and Ed) will 
simultaneously satisfy Eq. |A.2| for data at another frequency, for 
example, /y 2 , viz. 


A/v20-) = It - 2^d0-, j) Kriii) Bylina, ;•)], (A.5) 

where we let A/y 2 ( 7 ) ^ 0. Eq. |A.5| yields for each opacity curve, 
designated by 7 , a unique solution for T^{j) and Ed( 7 ). It only 
remains to select the correct /c-curve. Minimization of the data 
for a third frequency, /ys, could in principle be used to achieve 
this. Hence 


= Ky^U) - y / ( W • 

2d0)5y3[7^d0)] 

Minimizing A/Cy then provides the unique solution for the 
temperature Td, the surface density Ed and the opacity law of the 
dust, K - /c(y). Einally, the assumption of low optical depth at the 
different frequencies of this solution is then verified a posteriori 

(Eq.|g. 

The column density of the gas, A^g = N(H), needs to be de¬ 
termined from other observations. This, finally, leads to the gas- 
to-dust mass ratio from 


g2d = 


pmuN^ 

Sd 


(A.7) 


where p = 2.4 is the mean molecular weight. The derived dust 
opacity curve has been computed for a particular volume den¬ 
sity of the gas, ^^(H) and this ought to be consistent with that 
determined from the spectral line fitting, The spatial av¬ 

eraging is done over each pixel, so that the total dust mass, Md, 
is then obtained from the summation of the surface density per 
pixel over all pixels. 





